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1. Introduction 

This talk was given as a progress report summarizing our work on the perturbative gluon prop- 
agator that has, in the meantime, published in Ref. [1]. This work is a continuation of Ref. [2] that 
was focused on the ghost propagator. The Monte Carlo study of both propagators (for a recent dis- 
cussion see [3]), which are closely related to each other by Schwinger-Dyson equations (SDE) [4], 
has attracted much attention outside the lattice community by phenomenologists working in in- 
frared QCD and hadron physics. For many reasons, the (covariant) Landau gauge propagators have 
been the centre of interest. The importance for the confinement problem has been discussed in 
Ref. [5]. 

Taken together, both propagators provide us with a definition and the momentum dependence 
of the running coupling a s {q 2 ) directly based on the ghost-gluon vertex. This has recently turned 
out to be an interesting frame to accurately measure CC S (M Z ) from propagators [6, 7, 8]. 

A simple connection between the two propagators exists in the extreme infrared, powerlike in 
a scaling and massive in a decoupling version. Only the latter is found on the lattice, while both 
solutions can be accommodated in the SDE approach. Some ideas exist today [9] about how to 
handle the Gribov ambiguity such that the scaling solution could eventually be reproduced on the 
lattice. 

In a wider sense, the effect of nontrivial vacuum structure (vortices, instantons) is mani- 
fest [10, 11] also in the gluon propagator, in the intermediate momentum range around @{\ GeV) 
where the SDE approach suffers from truncation ambiguities and where nonperturbative lattice cal- 
culations are unrivalled. In order to understand the onset of nonperturbative effects, it is important 
to approach this momentum range from high momenta within higher-order perturbation theory. 
Some of us have started such a program a couple of years ago [12]. While ordinary diagrammatic 
lattice perturbation theory soon gets too involved to be pursued, Numerical Stochastic Perturbation 
Theory (NSPT, for a recent review see Ref. [13] and references therein), provides a powerful tool 
to perform high-loop computations. One has to run coupled Langevin simulations on the lattice 
and to perform the necessary limits: Langevin time step e — > 0, volume V — > °° and lattice spacing 

The infinite-volume and continuum extrapolation part of the program has been satisfactorily 
achieved in Refs. [2, 1], such that the otherwise difficult to access non-logarithmic contributions 
became calculable. We refer the interested reader to these papers for more details about the tech- 
nique of NSPT and the procedure to take the needed limits. The method to perform the V — > °° and 
a — > limits simultaneously has been outlined for the first time in Ref. [2] for the ghost propagator 
and applied to the gluon propagator in Ref. [1]. There we have attempted to compare the pertur- 
bative results summed to the presently known order with the results of Monte Carlo simulations. 
In this short account of our work we concentrate on the method to extract all non-leading-log and 
non-logarithmic coefficients of the gluon dressing function from NSPT data. 

2. The gluon dressing function 

Recalling that the gluon propagator has to be color-diagonal and symmetric in the Lorentz 
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indices, its tensor structure in the continuum can be written as 



Df v (p) = S ab 



(2.1) 



where the transverse part D(p 2 ) and the longitudinal part F(p 2 ) have been introduced. The latter 
vanishes in the Landau gauge. The lattice counterpart D c ^ v {p{k)) 1 of Eq. (2.1) generally contains 
more terms on the r.h.s. due to the lower symmetry. The quantity under study, surviving the 
continuum limit, is 

D(p(k))= l -^D^(p(k)), (2.2) 
- > Ai=l 



where color indices are dropped from now on to ease the notation. Eq. (2.2) obviously corresponds 



to the transverse propagator. In the following we use the w^-order dressing function Ji n \ (p) defined 



as 



J{n){p)=P 2 D (n) {p{k)) . (2.3) 

3. The non-logarithmic contributions to Z 

In the RF-MOM scheme, the relation between the bare dressing function Jo(p,a, OCgji) and its 
renormalized counterpart Jrp (/?,/!, a#/' ) is given by 2 

Jo(p,a,agj>) = Z(a,fi,a RI >)Jgp{p,fi,a ia i) , (3.1) 

with the renormalization condition Jrp (p, pi , 0Cri')\ p 2 = ^2 = 1 and with ot/y/ = g 2 {api) / \6n 2 . More 
in detail, the objects in Eq. (3.1) can be perturbatively expanded as 

n 

Z(a,n,asp) = 1 + £ cfo £ z^'j logman) , (3.2) 

n>0 j=0 

n / y \J 

J (p,a, aw) = 1 + £ a n m , £ ~ l°gM 2 , (3-3) 

„>o j=o V z / 

Jsrip^agr) = l + Z^'i^'jih^iP^ 1 ) 2 ) 1 ■ (3-4) 

„>o j=o V z / 

By plugging the perturbative expansions above into Eq. (3.1) and by converting the renormalized 
coupling into the bare one, Ofo = N c / (87T 2 j3), the bare dressing function can be written as 

Jo(p,a,p) = 1 + £ ^ £ J n j logi (ap) 2 , (3.5) 

n>0 P j=Q 



p(k) components are denned as Pn(kfi) = Ink^jL^, with = aNjx the lattice size along direction )l and k^ 6 
[0,Np — 1], In what follows, the subscript (ionA' and L will be dropped since simulations were performed on symmetric 
lattices. 

2 Formulae in this section closely mimic those in section 4 of [2] which we suggest the interested reader to consult. 
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where the coefficients J n j are related to the z„ ,'s in Eq. (3.2). As far as the logarithmic contribu- 
tions are concerned, their coefficients depend on the anomalous dimension of the gluon field and 
the j3 — function (see [14]). The purpose of this work is to compute the J n fiS which are related to 
thezjj's by 

7i,o = 0.03799544zfo , (3-6) 
J 2 ,o = 0.10673710zfJ + 0.00144365 zfj , (3.7) 
/ 3 ,o = 0.375990zfo + 0.008 11 105 zf 7 o + 0.0000548523 zf j . (3.8) 

4. The fitting procedure 

We can isolate the n th — loop contribution in Eq. (3.5) (without power of j3~") and write 

n 

J{n){p,a) = J n ,o(ap) + £ J n j log J (ap) 2 , (4.1) 
;'=i 

where, by recalling the existence of irrelevant lattice artifacts, J n o(ap) can be decomposed as 

(an) 4 

Jnfi{ap)=J nfl + c n . } {ap) 2 + c n ^j J ^+c n ^{apf + ... \(ap) m = £(t?/^) m ] , (4.2) 

\ a P) n 

where J n fi is the n — loop constant we want to compute. Taking into account also finite-size effects, 
J n fl(ap) has to be replaced by J n o(ap,pL) with 

J n ,o(ap,pL) = J n ,o(ap) + [J nfl (ap,pL) - J nfi (ap)} = J nfi (ap) + 8J nfl (ap,pL) = 

= J„fi(ap) + 8J n .o(0,pL) = J„. (ap) + 8J nfi (pL) , (4.3) 

where, in the last two steps, ap corrections on pL effects are assumed to be corrections on correc- 
tions and have been neglected. Furthermore we assume that the point p max corresponding to the 
largest (ap) 2 on the largest lattice size L max is such that 8J n) o(p maK L maK ) = 0. 

When treating pL contributions, it is useful to notice that, since p^L = p^aN = 2nk^, any 
given 4-tuple of integers (&i,&2)&3)&4) nas the same finite-size corrections to J n ,o(ap) on different 
lattice sizes. Therefore, the steps in data-analysis at any loop order n can be summarized as follows: 



• a suitable window [(a/?)^ in , (ap)m ax ] is identified where a sufficient number of data points is 
available on different lattice sizes; 

• the logarithmic contributions are subtracted from the bare dressing function in order to get 

J„ } o(ap,pL); 

• J n .o(ap,pL) is fitted according to the procedure sketched above: from Eq. (4.2) we compute 
the desired J Ki q. 



Obviously, only results stemming from stable fits have to be taken into account. 
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5. Results 



Figs. 1, 3 and 5 show an example of the fitting procedure outlined above: the black dots 
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log subtracted data ■ 
no pL 
no pL, H{4) 
final fit 2.298(40) 



1-loop 

X 2 /d.o.f. =0.9562 
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2.45 - 
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2.25 - 
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2.15 - 



exact 2.29368 
(J lfi ) = 2.303(34) 



(apY 
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Figure 1: 1-loop data vs. (ap) 2 : the different sets Figure 2: J\ o results vs. x 2 '- the dotted blue line 
of points stem from the various steps of the analy- is the average of the 10 blue points, the continuous 



sis (see section 4). 



in 



9 - 



log subtracted data 
no pL 
no pL. H(4) 
final fit 8.032(122) 

2-loop 

X 2 /d.o.f.=1.0633 



i i It 1 2 



red one is the analytical result. 

' (J 2 ,„>' = 7.921(117) 



4 6 

(ap) 2 



Figure 3: Same as in Fig. 1 for 2-loop data. 
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Figure 4: J2,o results vs. % 2 : see caption of Fig. 2. 



stand for log-subtracted data, blue stars for data after finite-size effects have been removed while 
red spots for data after further removal of hypercubic effects. 



Loop order n 


1 


2 


3 


Jn,0 


2.30(3) 


7.92(12) 


31.7(5) 



Table 1: Results for the constant contributions to Z up to the third loop. Compare the 1-loop analytical 
result h o = 2.29368. 



In general, there could be various fitting functions that fulfill the stability requirement men- 
tioned at the end of section 4: the results we quote come from averaging the values obtained from 
the fits with the 10 lowest % 2 at every perturbative order: Figs. 2, 4 and 6 contain these values 
and their average with their errorbars. The final outcome of the analysis sketched in the above is 
contained in Table 1 . 
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log subtracted data 
no pL 
no pL. HU) 
final fit 31.53(56) 

3-loop 

x 7d.oi.=1.1085 



12 3 4 



(J 3>0 ) = 31.66(47) 



6 7 



(ap) 2 

Figure 5: Same as in Fig. 1 for 3-loop data. 
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Figure 6: ^3,0 results vs. % 2 : see caption of Fig. 2. 



Resumming the perturbative series for Jo(a,p,fi) up to the third loop, we get 



J 3 - loop (p,a,j5) = 1+^ (-0.24697 log(a/?) 2 + 2.29368^ + 
H 

1 



+ — (0.08211 (\og(ap) 2 ) 2 - 1. 48445 log(a^) 2 + 7.93(12)) + (5.1) 



+ 



0.02964 (log(ap) 2 ) +0.81689 (\og{ap) 2 ) + 



1 

w 

8.13(3) \og{ap) 2 + 31.7(5) 



which can then be converted to the RF-MOM scheme by using the formulae of section 3. 
6. Conclusion 

In this work we presented an algorithm to compute the non-logarithmic parts of infrared di- 
vergent quantities in NSPT in the infinite- volume limit using the example of the gluon propagator. 
The one-loop result coincides with the well-known result from the diagrammatic approach [15]. 
The corresponding two- and three-loop finite constants have been computed for the first time. 

We can use the dressing function obtained with NSPT at finite lattice sizes to investigate the 
perturbative background of the corresponding Monte-Carlo calculated propagators. This has been 
done in some detail in [1]. 
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